// A digital media literacy intervention increases 
// discernment between mainstream and false news in the United States and India)
// Guess et al.
//
// Replication Code for India Results (Replication Archive)

*Make sure cibar, estout, lincomest, lean1, ivreg2, ranktest, balancetable are installed

clear
cd "\\afs\umich.edu\user\l\e\lernerm\Cleaned India Data and Code" // Set working directory

//// Appendix Balance
// Table B1
// Perceived accuracy of false headlines for which fact-check messages were distributed (Wave 2)
**** Online
use "Data/onlineRQE1w2-onlineRQE2w2_data.dta", clear
reg accuracyw2 i.control_placebo_factcheck##tips i.dv if false==1 & factchecked==1 & control_placebo_factcheck !=0, robust cluster(caseid)
est store onlineRQA1w2
lincomest 2.control_placebo_factcheck+2.control_placebo_factcheck#1.tips
est store onlineRQA1w2_lincom

**** Face-to-face
use "Data/ftfRQE1w2-ftfRQE2w2_data.dta", clear
reg accuracyw2 i.placebo_factcheck##tips i.dv if false==1 & factchecked==1, robust cluster(caseid)
est store ftfRQA1w2
lincom 2.placebo_factcheck+2.placebo_factcheck#1.tips
est store ftfRQA1w2_lincom

** Output Table B1
estout onlineRQA1w2 ftfRQA1w2, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 
estout onlineRQA1w2_lincom ftfRQA1w2_lincom, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 

// Table B2
// Effects of India fact-check message experiment on perceived headline accuracy
**** Online
use "Data/onlineRQE1w2-onlineRQE2w2_data.dta", clear
reg accuracyw2 i.control_placebo_factcheck i.dv if false==1 & factchecked==1 & control_placebo_factcheck != 0, robust cluster(caseid)
est store onlineHB1aw2
reg accuracyw2 i.control_placebo_factcheck##factchecked if false==1 & control_placebo_factcheck != 0, robust cluster(caseid)
est store onlineRQB1w2

**** Face-to-face
use "Data/ftfRQE1w2-ftfRQE2w2_data.dta", clear
reg accuracyw2 i.placebo_factcheck i.dv if false==1 & factchecked==1, robust cluster(caseid)
est store ftfHB1aw2
reg accuracyw2 i.placebo_factcheck##factchecked if false==1, robust cluster(caseid)
est store ftfRQB1w2

** Output Table B2
estout onlineHB1aw2 onlineRQB1w2 ftfHB1aw2 ftfRQB1w2, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 

// Table B3
// Effect of India fact-check message experiment on perceived accuracy of false news head-line accuracy for which fact-check messages were distributed (Wave 2) by partisan congeniality
**** Online
use "Data/onlineRQE1w2-onlineRQE2w2_data.dta", clear
reg accuracyw2 i.control_placebo_factcheck##congenial i.dv if false==1 & factchecked==1 & control_placebo_factcheck != 0, robust cluster(caseid)
est store onlineHC1bw2

**** Face-to-face
use "Data/ftfRQE1w2-ftfRQE2w2_data.dta", clear
reg accuracyw2 i.placebo_factcheck##congenial i.dv if false==1 & factchecked==1, robust cluster(caseid)
est store ftfHC1bw2

** Output Table B3
estout onlineHC1bw2 ftfHC1bw2, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)

// Table B4
// Effect of India fact-check experiment on perceived accuracy of false fact-checked newsheadline accuracy by partisan congeniality (Wave 2, question-level)
**** Online
use "Data/onlineRQE1w2-onlineRQE2w2_data.dta", clear
levelsof dv if false==1 & factchecked==1, local(dvs)
foreach i of local dvs {
     reg accuracyw2 i.control_placebo_factcheck##bjp_support if false==1 & factchecked==1 & control_placebo_factcheck != 0 & dv == `i', robust cluster(caseid)
	est store online_subHC1bw2_`i'
}

**** Face-to-face
use "Data/ftfRQE1w2-ftfRQE2w2_data.dta", clear
levelsof dv if false==1 & factchecked==1, local(dvs)
foreach i of local dvs {
     reg accuracyw2 i.placebo_factcheck##bjp_support if false==1 & (bjp_support==1|bjp_oppose==1) & factchecked==1 & dv==`i', robust cluster(caseid)
	est store ftf_subHC1bw2_`i'
}
** Output Table B4
estout online_subHC1bw2_* ftf_subHC1bw2_*, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)

// Table B5
// Effect of india media literacy intervention on perceived accuracy of false headlines bysample and fact-checking condition
use "Data/pooledthreeway_data.dta", clear
*ITT
reg accuracyw2 tips##factcheck_condition##ftf i.dv, robust cluster(caseid)
est store pooledthreeway

*ATT
ivreg2 accuracyw2 (r_tips r_tips_x_fc_condition r_tips_x_ftf r_tips_x_fc_condition_x_ftf = tips tips_x_fc_condition tips_x_ftf tips_x_fc_condition_x_ftf) ftf factcheck_condition fc_condition_x_ftf i.dv, robust cluster(caseid)
est store pooledthreeway_att

** Output Table B5
estout pooledthreeway pooledthreeway_att, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 

// Table B6
// Association of partisan congeniality on perceived accuracy
**** Online
use "Data/onlineHA1-onlineHC1_data.dta", clear
reg accuracy bjp_support##pro_bjp##anti_bjp bjp_oppose##pro_bjp##anti_bjp i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest if false==1, robust cluster(caseid)
est store onlineHD1
lincomest 1.bjp_support#1.pro_bjp-1.bjp_oppose#1.pro_bjp
est store onlineHD1_lincom_probjp
reg accuracy bjp_support##pro_bjp##anti_bjp bjp_oppose##pro_bjp##anti_bjp i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest if false==1, robust cluster(caseid)
lincomest 1.bjp_support#1.anti_bjp-1.bjp_oppose#1.anti_bjp
est store onlineHD1_lincom_antibjp

**** Face-to-face
use "Data/ftfHA1-ftfHC1_data.dta", clear
reg accuracy bjp_support##pro_bjp##anti_bjp bjp_oppose##pro_bjp##anti_bjp i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest if false==1, robust cluster(caseid)
est store ftfHD1
lincomest 1.bjp_support#1.pro_bjp - 1.bjp_oppose#1.pro_bjp
est store ftfHD1_lincom_probjp
reg accuracy bjp_support##pro_bjp##anti_bjp bjp_oppose##pro_bjp##anti_bjp i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest if false==1, robust cluster(caseid)
lincomest 1.bjp_support#1.anti_bjp - 1.bjp_oppose#1.anti_bjp
est store ftfHD1_lincom_antibjp

** Output Table B6
estout onlineHD1 ftfHD1, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 
estout onlineHD1_lincom_probjp ftfHD1_lincom_probjp, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 
estout onlineHD1_lincom_antibjp ftfHD1_lincom_antibjp, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 

// Table B7
// Association of nationalism and perceived accuracy of nationalistic false news headlines
**** Online
use "Data/onlineHA1-onlineHC1_data.dta", clear
reg accuracy indian_defense_goal##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store onlineRQD11
reg accuracy c.indian_use_military##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store onlineRQD12
reg accuracy c.indian_pakistan_threat##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store onlineRQD13
reg accuracy c.hindu_ram_mandir##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store onlineRQD14
reg accuracy c.hindu_muslim_pop##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store onlineRQD15

**** Face-to-face
use "Data/ftfHA1-ftfHC1_data.dta", clear
reg accuracy indian_defense_goal##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store ftfRQD11
reg accuracy c.indian_use_military##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store ftfRQD12
reg accuracy c.indian_pakistan_threat##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store ftfRQD13
reg accuracy c.hindu_ram_mandir##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store ftfRQD14
reg accuracy c.hindu_muslim_pop##nationalist i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose if false==1, robust cluster(caseid)
est store ftfRQD15

** Output Table B7
estout onlineRQD1*, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 
estout ftfRQD1*, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 

// Table B8
// Association of media trust and perceived accuracy of false news headlines
**** Online
use "Data/onlineHA1-onlineHC1_data.dta", clear
reg accuracy c.newspaperaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store onlineRQD21
reg accuracy c.televisionaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store onlineRQD22
reg accuracy c.facebookaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store onlineRQD23
reg accuracy c.whatsappaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store onlineRQD24

**** Face-to-face
use "Data/ftfHA1-ftfHC1_data.dta", clear
reg accuracy c.newspaperaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store ftfRQD21
reg accuracy c.televisionaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store ftfRQD22
reg accuracy c.facebookaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store ftfRQD23
reg accuracy c.whatsappaccuracy i.agegroup i.caste i.education i.female i.religion political_knowledge political_interest bjp_support bjp_oppose i.dv if false==1, robust cluster(caseid)
est store ftfRQD24

** Output Table B8
estout onlineRQD2* ftfRQD2*,replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex) 

//Table B9
// Effect of India media literacy intervention on perceived accuracy by partisan congeniality (question-level)
**** Online
use "Data/onlineHA1-onlineHC1_data.dta", clear
levelsof dv if false==1 , local(dvs)
foreach i of local dvs {
     reg accuracy tips##bjp_support if false==1 & (bjp_support==1|bjp_oppose==1) & dv==`i', robust cluster(caseid)
	est store online_subHC1a1_`i'
}

levelsof dv if true==1 , local(dvs)
foreach i of local dvs {
     reg accuracy tips##bjp_support if true==1 & (bjp_support==1|bjp_oppose==1) & dv==`i', robust cluster(caseid)
	est store online_subHC1a2_`i'
}

**** Face-to-face
use "Data/ftfHA1-ftfHC1_data.dta", clear
levelsof dv if false==1, local(dvs)
foreach i of local dvs {
    reg accuracy tips##bjp_support if false==1 & (bjp_support==1|bjp_oppose==1) & dv == `i', robust cluster(caseid)
	est store ftf_subHC1a1_`i'
}

levelsof dv if true==1, local(dvs)
foreach i of local dvs {
    reg accuracy tips##bjp_support if true==1 & (bjp_support==1|bjp_oppose==1) & dv == `i' , robust cluster(caseid)
	est store ftf_subHC1a2_`i'
}

** Output Table B9
estout online_subHC1a1_* ftf_subHC1a1_*, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)
estout online_subHC1a2_* ftf_subHC1a2_*, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)


// Table B10
// Effect of India media literacy intervention on perceived accuracy by Indian and Hindu nationalism
**** Online
use "Data/onlineHA1-onlineHC1_data.dta", clear
reg accuracy tips##indian_defense_goal i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_1
reg accuracy tips##c.indian_use_military i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_2
reg accuracy tips##c.indian_pakistan_threat i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_3
reg accuracy tips##indian_defense_goal i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_4
reg accuracy tips##c.indian_use_military i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_5
reg accuracy tips##c.indian_pakistan_threat i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_6
reg accuracy tips##c.hindu_ram_mandir i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_7
reg accuracy tips##c.hindu_muslim_pop i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_8
reg accuracy tips##c.hindu_ram_mandir i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_9
reg accuracy tips##c.hindu_muslim_pop i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store onlineRQC1_10

**** Face-to-face
use "Data/ftfHA1-ftfHC1_data.dta", clear
reg accuracy tips##indian_defense_goal i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_1
reg accuracy tips##c.indian_use_military i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_2
reg accuracy tips##c.indian_pakistan_threat i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_3
reg accuracy tips##indian_defense_goal i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_4
reg accuracy tips##c.indian_use_military i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_5
reg accuracy tips##c.indian_pakistan_threat i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_6
reg accuracy tips##c.hindu_ram_mandir i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_7
reg accuracy tips##c.hindu_muslim_pop i.dv if false==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_8
reg accuracy tips##c.hindu_ram_mandir i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_9
reg accuracy tips##c.hindu_muslim_pop i.dv if true==1 & nationalistic==1, robust cluster(caseid)
est store ftfRQC1_10

** Output Table B10
estout onlineRQC1_1 onlineRQC1_4 onlineRQC1_2 onlineRQC1_5 onlineRQC1_3 onlineRQC1_6 onlineRQC1_7 onlineRQC1_9 onlineRQC1_8 onlineRQC1_10, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)
estout ftfRQC1_1 ftfRQC1_4 ftfRQC1_2 ftfRQC1_5 ftfRQC1_3 ftfRQC1_6 ftfRQC1_7 ftfRQC1_9 ftfRQC1_8 ftfRQC1_10, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)

// Table B11
// Effect of India media literacy intervention on perceived accuracy by Indian and Hindu nationalism (question-level)
**** Online
use "Data/onlineHA1-onlineHC1_data.dta", clear
levelsof dv if false==1 & nationalistic==1 , local(dvs)
foreach i of local dvs {
	reg accuracy tips##indian_defense_goal if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_1_F_`i'
	reg accuracy tips##c.indian_use_military if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_2_F_`i'
	reg accuracy tips##c.indian_pakistan_threat if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_3_F_`i'
	reg accuracy tips##c.hindu_ram_mandir if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_4_F_`i'
	reg accuracy tips##c.hindu_muslim_pop if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_5_F_`i'
}

levelsof dv if true==1 & nationalistic==1, local(dvs)
foreach i of local dvs {
	reg accuracy tips##indian_defense_goal if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_1_T_`i'
	reg accuracy tips##c.indian_use_military if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_2_T_`i'
	reg accuracy tips##c.indian_pakistan_threat if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_3_T_`i'
	reg accuracy tips##c.hindu_ram_mandir if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_4_T_`i'
	reg accuracy tips##c.hindu_muslim_pop if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store online_subRQC1_5_T_`i'
}

**** Face-to-face
use "Data/ftfHA1-ftfHC1_data.dta", clear
levelsof dv if false==1 & nationalistic==1 , local(dvs)
foreach i of local dvs {
	reg accuracy tips##indian_defense_goal if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_1_F_`i'
	reg accuracy tips##c.indian_use_military if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_2_F_`i'
	reg accuracy tips##c.indian_pakistan_threat if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_3_F_`i'
	reg accuracy tips##c.hindu_ram_mandir if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_4_F_`i'
	reg accuracy tips##c.hindu_muslim_pop if false==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_5_F_`i'
}

levelsof dv if true==1 & nationalistic==1, local(dvs)
foreach i of local dvs {
	reg accuracy tips##indian_defense_goal if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_1_T_`i'
	reg accuracy tips##c.indian_use_military if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_2_T_`i'
	reg accuracy tips##c.indian_pakistan_threat if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_3_T_`i'
	reg accuracy tips##c.hindu_ram_mandir if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_4_T_`i'
	reg accuracy tips##c.hindu_muslim_pop if true==1 & nationalistic==1 & dv == `i', robust cluster(caseid)
	est store ftf_subRQC1_5_T_`i'
}

** Output Table B11
estout online_subRQC1_1_F_5 online_subRQC1_1_F_6 online_subRQC1_1_T_11 online_subRQC1_1_T_12 online_subRQC1_2_F_5 online_subRQC1_2_F_6 online_subRQC1_2_T_11 online_subRQC1_2_T_12 online_subRQC1_3_F_5 online_subRQC1_3_F_6 online_subRQC1_3_T_11 online_subRQC1_3_T_12, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)
estout online_subRQC1_4_F_5 online_subRQC1_4_F_6 online_subRQC1_4_T_11 online_subRQC1_4_T_12 online_subRQC1_5_F_5 online_subRQC1_5_F_6 online_subRQC1_5_T_11 online_subRQC1_5_T_12, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)
estout ftf_subRQC1_1_F_5 ftf_subRQC1_1_F_6 ftf_subRQC1_1_T_11 ftf_subRQC1_1_T_12 ftf_subRQC1_2_F_5 ftf_subRQC1_2_F_6 ftf_subRQC1_2_T_11 ftf_subRQC1_2_T_12 ftf_subRQC1_3_F_5 ftf_subRQC1_3_F_6 ftf_subRQC1_3_T_11 ftf_subRQC1_3_T_12, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)
estout ftf_subRQC1_4_F_5 ftf_subRQC1_4_F_6 ftf_subRQC1_4_T_11 ftf_subRQC1_4_T_12 ftf_subRQC1_5_F_5 ftf_subRQC1_5_F_6 ftf_subRQC1_5_T_11 ftf_subRQC1_5_T_12, replace varwidth(25) collabels("") cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) stats(r2 N, fmt(%9.2f %9.0f) labels("R^2" "N")) starlevels(* 0.05 ** 0.01 *** 0.005) style(tex)